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Here we respond to the comment by Tsagas [ij on our paper J^] . We show that the results in that 
comment are flawed and cannot be used for drawing conclusion about the nature of magnetic field 
amplification by gravitational waves, and give further support that the results of 01 are correct. 
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I. INTRODUCTION 

In Ref. ^ it was shown, using proper second order 
covariant gauge-invariant perturbation theory, that the 
results in Ref. ^] concerning gravitational wave ampli- 
fication of weak magnetic fields gave amplification rates 
incorrect by several orders of magnitude. The author 
of the comment pL] claims that our main result in Ref. 
does not apply outside the Hubble radius, and that 
our numerical estimates and conclusions are therefore in- 
valid as a result of neglecting oscillatory terms. Further 
claims that our approach is not truly gauge-invariant and 
also not mathematically complete as a result of neglect- 
ing the issue of constraints are also made. However, all 
these claims will be refuted in what follows. Moreover, 
the comment contains misleading statements (address- 
ing some of our assumptions and approximations) which 
are in contrast to what was done (and indeed clearly 
stated) in our paper. In particular, we did not employ 
ideal magnetohydrodynamics as the comment states; it 
was only assumed that the electric field at first order 
is perturbatively smaller than the magnetic field. That 
the gravito-magnetic interaction did not yield an electric 
field at second order is a consequence of requiring an ho- 
mogeneous magnetic field interacting with gravity waves. 
Ideal magnetohydrodynamics was only used in Ref. 



II. THE ISSUE OF THE ELECTRIC CURL 

In section II of the comment [l| , the author argues that 
because of Eq. (1) it is not possible a priori to assume 
that (curl Ea)^ = initially even though curl Ea =0 ini- 
tially. However, it follows from Maxwell's equations that 
(curli?a)'_L = —Qz\ir\ Ea -f curl curl and since initially 
curlSa vanishes (and therefore also curl curl S^,) because 



the first-order magnetic field Ba is homogeneous, one 
concludes that also (curlE'^)^ = when the interaction 
between gravitational waves and the first-order magnetic 
field is turned on. When these constraints are used in 
the governing equation for curli^a (see Eq. 2 in Q), one 
indeed finds that in the spatially flat case the generated 
electric field stays curl-free if it was initially so. 

The author of the comment came to the opposite con- 
clusion by using the relation 



curl curl B„ 



(1) 



However, Ba is not a gauge-invariant variable at second 
order, and thus cannot be used for the purpose of analyz- 
ing the appropriateness of the initial conditions. Let us 
make this important point very clear: the usual splitting 
of the magnetic field Ba into a first-order homogeneous 

(2) 

part and a second-order term, Ba = Ba + Ba , as the 
author of the comment assumed, leads to critical incon- 
sistencies when commutation relations are involved (see 
also Eqs. (15) and (16) in 0). Taking the spacetime 
background to be spatially flat and inserting the said 
splitting in the right hand side of Eq. ([1]) yields 



curl curl Ba = -D^B^^) 



(2) 



where TZ^^J denotes the first-order contribution to the 3- 
Ricci tensor. On the other hand, if one first employs the 
splitting for the term curl Ba, and then applies Eq. ([1]), 
one obtains 



curl curl Ba = curl curl Bi^-* 



-D'B(p, (3) 



*Also at: Centre for Fundamental Physics, Rutherford Applcton 
Laboratory, Chilton, Didcot, Oxon OXll OQX, U.K. 



which differs from the result This ambiguity, due 
to the gauge problem with the magnetic field splitting in 
relation to the use of the commutator relation, renders 
Eq. (1) in [1] meaningless and invalidates the conclusion 
drawn from it. 

We stress that we neglected the electric current at all 
orders in our paper , which is potentially physically un- 
sound. However, the inclusion of a first-order current re- 
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quires the inclusion of a first-order inhomogeneous mag- 
netic field, which was beyond the scope of the paper Q . 
Currents and velocity perturbations have been explicitly 
taken into account in our follow-up paper [sj by making 
use of ideal magnctohydrodynamics. In this case we did 
indeed have {curlE"") ^ 0, it was shown that this term 
only contributes at very small scales to the generated 
magnetic field. On large scales, the results for the gener- 
ated magnetic field agree with the ones found previously 
(compare Eq. (50) of [s^ with Eq. (49) of [2]). 



III. SCALES AND AMPLIFICATION 
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Dust: L=0.01 
Radiation: L=0.01 



When discussing super-horizon scales in section III 
of his comment, the author correctly points out that 
the solution Eq. (7) in Ref. [2] of the interaction 
term J(t) in the case of dust is valid as long as a; = 
2£r^/^/(aoiJo) ^ 1 holds (similarly for the case of ra- 
diation). However, since the gravitational wave num- 
ber is defined a& i = 2TTa/XGw and the Hubble length 
\h = ^/H, we also have x — 2Tr{XH/XGw)oT^^^, and 
therefore the condition a; ^ 1 will sooner or later break 
down as r grows, even though the ratio {Xh / ^gw)o is 
small for a particular gravity wave mode with wavenum- 
ber £. The solutions Eq. (4) and Eq. (5) in [H therefore 
describe the evolution of an initially super-horizon mag- 
netic mode [meaning {Xh / Xgw)o ^ 1] only correctly up 
to the time of horizon crossing (at best). Thereafter one 
has to rely on the solutions Eqs. (49) and (50) given in 

a. 

The fact that we did not display the solutions of the 
generated magnetic field for the dust and radiation case 
in full detail when dealing with the modes with non-zero 
wavenumber, £ ^ 0, seems to have led to some confu- 
sion. To arrive at the solutions Eqs. (49) and (50) in 
we used generic initial conditions leading to complicated 
expressions, and they were thus not stated explicitly in 
Q; however, all the relevant terms for the amplification 
have been given and the nature of the non-displayed ex- 
pressions has been circumscribed in 0]. To improve on 
this, we give in the following the full solution for the total 
(first- plus second-order) magnetic field in the radiation 

era assuming initial conditions I(i)(t = 1) = a^^^ Bq and 
/^'^^(l) = for the interaction variable. The latter con- 
dition simply means that ((T(f)/iJ)o = 2(Tq^'' has been 
chosen. With these initial conditions, Eq. (50) in Q 
gives 
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FIG. 1: The dimensionless induced magnetic fields plotted 
against dimensionless time r for dust and radiation, taking 
L = 0.01. For given initial conditions the generated mag- 
netic field saturates at a value described by our result ([6}. 
The super-adiabatic growth is approximately depicted by the 
infinite- wavelength solutions (46) and (47) in Q]. 



where L = £/{aoHo) — 2'k{Xh / Xgw)o- As expected, 
the generated magnetic field is proportional to the ini- 
tial shear anisotropy {a / H){)] it starts to grow from a 
zero value until it saturates to a constant value (mod- 
ulo adiabatic decay); this is also evident from FIG. [T] 
Note that the solution (jl]) is exact and valid for all finite 
wavelengths at all times. In the infinite-wavelength limit, 
£ — > 0, the solution for the magnetic field goes over into 



<a(«) _ 




- 1 



(5) 



which is of course the same as Eq. (47) in 0, where 
identical initial conditions had been used. 

A numerical analysis of the gravito-magnetic interac- 
tion (cf. FIG. H]) demonstrates that the generated mag- 
netic field starts from zero, undergoes super-adiabatic 
growth until it begins to oscillate around a value which 
depends on the initial conditions but is the same dur- 
ing both the radiation or dust/reheating eras. Whilst 
the super-adiabatic growth is approximately described 
by the infinite-wavelength solutions Eqs. (46) and (47) 
in 0, the saturation effect is only present in the finite- 
wavelength solutions Eqs. (49) and (50) in [2] [see also 
Eq. (g]) above]. 

What then is the emerging magnetic field at the end 
of the day? Returning to Eq. (g]) and dividing by the 
energy density of the background radiation, one obtains 
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(the oscillatory parts decay away) 
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(6) 



where the wavenumber indices have been suppressed and 
the resonant condition Aqw ^ — gravity wave- 
length Agw matches the size A^ of the magnetic region 
— has been used. This agrees with our result Eq. (51) 
given in 0. As pointed out in our paper, the case of dust 
is very similar and leads to the same result ([6]) for the 
maximally resulting magnetic field. 

It is plain to see that in the light of the above (being in 
unison with comments already made in our paper 0), the 
criticism brought forward by the author of the comment 
(see Eqs. (8) and (9) therein) is not correct. Also, Eq. 

(9) of [1] is algebraically wrong and can not be used to 
draw conclusions about the validity of the results in . 

In the last paragraph of section III of the comment, the 
author claims that our numerical results have been com- 
promised because we did not take reheating effects into 
account. He further claims that one cannot use the result 
Eq. ([6]) for this purpose and presents instead formula Eq. 

(10) in [ij]. However, that formula is of third order, and 
is therefore not correct since the gravito-magnetic inter- 
action is of second order. If one assumes that reheating 
does not last long enough in order for the induced mag- 
netic field to reach saturation, one may use Eq. (5) of [1] 
to find the field at the end of reheating, which may then 
be employed in our result Eq. ([6]) to yield the total mag- 
netic field during the radiation era. Including as always 
only terms up to second order, we obtain 
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correcting Eq. (10) of [IJ]. 

Contrary to the claim of the comment, one can very 
well use our result © , which was shown in to hold for 
the radiation as well as the dust /reheating era, in order to 
take reheating into account. Let us for simplicity assume 
that reheating lasts long enough such that the generated 
magnetic field saturates during the reheating phase. In a 
first step, we employ the result ^ to find the total field 
at the end of reheating, 



Bo \ aRH / 
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(8) 



Here, the suffix denotes the time at the end of infiation 
when reheating starts. In a second step we use again 
Eq. ^ but this time with Brh ^ input magnetic field. 
We then find that the total magnetic field during the 
radiation era is given by 



B{a) 
Bo 



1 

lb 



Ah 



1 

lb 



Ah 



where a third-order term has been discarded. In our 
article the second gravitational wave term in square 
brackets in © was carefully calculated and found to be 
of order ^ 10~^, as expected from a perturbation calcu- 
lation [5i] . The first gravitational wave term in the square 
bracket represents the reheating correction and since it 
is of the same perturbative order as the term previously 
discussed, its value will not be much different (in any case 
< 1). The same comment applies when one uses Eq. ([7]) 
instead of Eq. Q . The claim that our numerical results 
have been compromised is thus refuted. 



IV. GAUGE-INVARIANCE AND LINEARITY 

The author of the comment claims that obtaining the 
second-order generated magnetic field via integration of 
the gauge-invariant variable Pa — B^a> + 2QBa/3 is 
in itself not a gauge-invariant procedure, 
we expand the magnetic field Ba = B^^ 
into first, second and higher order parts. 



RH^i?/RH 

(9) 



tain Pa — B 



- r(2) 



However, if 
then we ob- 

2QBa /3, and /3a thus evidently only 
describes the second-order part (higher order terms are 

neglected as usual). To get Ba from an integration is 
harmless since, as we pointed out in section II above and 
also in section IIC of our paper , the gauge issue arises 
if commutator relations have to be calculated, which is 
not the case here. 

The author of the comment further writes in the sec- 
ond paragraph of section IV [l| that there was no new 
information in our second-order approach in comparison 
to his own first-order approach to the problem at hand 
This makes little sense, since the two methods are 
fundamentally different in nature, and also since the gov- 
erning equations for the generated magnetic field differ 
in both cases. Moreover, when solving numerically the 
relevant weak-field equations [see Eqs. (61) and (62) in 
2]] for the magnetic field, taking identical initial condi- 
tions as discussed above, one obtains different results (cf. 
FIG. [2]). In contrast to the magnetic field obtained in 
our second-order approach (see FIG. [1]), the weak-field 
magnetic field shows large oscillations but no saturation. 



V. SUMMARY 

In this response to the criticism brought forward by 
the comment T| of Tsagas, we pointed out the following: 

• The author's argumentation regarding the curl of the 
electric field is untenable since it is based on a relation 
(Eq. (1) in [Tj) which is not valid at second order; 

• The setting we employed was not that of ideal mag- 
netohydrodynamics ; 

• Further strengthening our main results in ^ with 
numerical experiments demonstrates that there was no 
improper assessment of scale in 0] ; 

• Equation (9) in l] is algebraically wrong, and can 
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FIG. 2: The generated magnetic field in the weak-field ap- 
proximation in dimensionless units. The induced magnetic 
fields are plotted against the dimensionless time variable r 
for the dust and radiation cases taking L = £/{aoHo) — 
27r(Ai//AGw)o = 0.01. 



thus not be used for drawing conclusions about the va- 
hdity of the expressions in [2] ; 

• Although our formalism explicitly makes use of sec- 
ond order perturbation theory, the author of the com- 
ment claims that our numerical results are comprised 
because we neglected a third-order term [see (10) in fij]; 

• The generated magnetic field obtained by means of 
the weak- field approximation (see Jj) differs significantly 
from the one obtained with our formalism, the former 
shown to give rise to unreasonable amplification rates; 

• The second-order variable f3a = B^a> + '2'dBa/3 con- 

(2) 

tains only the second-order magnetic field Ba , and that 
an integration of (3a with respect to time yields the evo- 
lution of B^a^ . 

In conclusion, the criticism contained in the com- 
ment has thus been shown to be unfounded. Rather, 
it should be pointed out that, if treated using the proper 
gauge-invariant covariant second order perturbation the- 
ory, the interaction between magnetic fields and gravita- 
tional waves can give rise to interesting effects and be a 
possible source of a boost of weak seed fields as energy is 
transfered between the different degrees of freedom. 
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